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A primary objective of the Lunar Laser Ranging (LLR) experiment is to provide precise 
observations of the lunar orbit that contribute to a wide range of science investigations. 
In particular, time series of the highly accurate measurements of the distance between the 
Earth and Moon provide unique information used to determine whether, in accordance 
with the Equivalence Principle (EP), both of these celestial bodies are falling towards the 
Sun at the same rate, despite their different masses, compositions, and gravitational self- 
energies. 35 years since their initiation, analyses of precision laser ranges to the Moon 
continue to provide increasingly stringent limits on any violation of the EP. Current 
LLR solutions give (—1.0 ± 1.4) X 10~ 13 for any possible inequality in the ratios of the 
gravitational and inertial masses for the Earth and Moon, A(Mq /Mj). This result, in 
combination with laboratory experiments on the weak equivalence principle, yields a 
strong equivalence principle (SEP) test of A(M g /M/)sep = (-2.0 ±2.0) X 10~ 13 . Such 
an accurate result allows other tests of gravitational theories. The result of the SEP test 
translates into a value for the corresponding SEP violation parameter r\ of (4.4 ± 4.5) X 
10 -4 , where r\ = 4/3 — 7 — 3 and both 7 and /3 are parametrized post-Newtonian (PPN) 
parameters. Using the recent result for the parameter 7 derived from the radiometric 
tracking data from the Cassini mission, the PPN parameter f) (quantifying the non- 
linearity of gravitational superposition) is determined to be /3 — 1 = (1.2 ± 1.1) X 10 — 4 . 
We also present the history of the lunar laser ranging effort and describe the technique 
that is being used. Focusing on the tests of the EP, we discuss the existing data, and 
characterize the modeling and data analysis techniques. The robustness of the LLR 
solutions is demonstrated with several different approaches that are presented in the text. 
We emphasize that near-term improvements in the LLR ranging accuracy will further 
advance the research of relativistic gravity in the solar system, and, most notably, will 
continue to provide highly accurate tests of the Equivalence Principle. 

Keywords: lunar laser ranging; equivalence principle; tests of general relativity. 



1. Introduction 

The Equivalence Principle (EP) has been a focus of gravitational research for more 
than four hundred years. Since the time of Galileo (1564-1642) it has been known 
that objects of different mass and composition accelerate at identical rates in the 
same gravitational field. In 1602-04 through his study of inclined planes and pen- 
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dulums, Galileo formulated a law of falling bodies that led to an early empirical 
version of the EP. However, these famous results would not be published for an- 
other 35 years. It took an additional fifty years before a theory of gravity that 
described these and other early gravitational experiments was published by Newton 
(1642-1727) in his Principia in 1687. Newton concluded on the basis of his second 
law that the gravitational force was proportional to the mass of the body on which 
it acted, and by the third law, that the gravitational force is proportional to the 
mass of its source. 

Newton was aware that the inertial mass Mj appearing in the second law 
F = M/a, might not be the same as the gravitational mass Mq relating force 
to gravitational field F = Mag. Indeed, after rearranging the two equations above 
we find a = [Mq / '.Mj)g and thus in principle materials with different values of the 
ratio (Mg/Mi) could accelerate at different rates in the same gravitational field. 
He went on testing this possibility with simple pendulums of the same length but 
different masses and compositions, but found no difference in their periods. On 
this basis Newton concluded that (Mg/Mi) was constant for all matter, and by a 
suitable choice of units the ratio could always be set to one, i.e. (Mg/Mi) = 1. 
Bessel (1784-1846) tested this ratio more accurately, and then in a definitive 1889 
experiment Eotvos was able to experimentally verify this equality of the inertial 
and gravitational masses to an accuracy of one part in 10 9 (see Refs. HTJ 142 } fl"2|) . 

Today, almost three hundred and twenty years after Newton proposed a compre- 
hensive approach to studying the relation between the two masses of a body, this 
relation still continues to be the subject of modern theoretical and experimental 
investigations. The question about the equality of inertial and passive gravitational 
masses arises in almost every theory of gravitation. Nearly one hundred years ago, 
in 1915, the EP became a part of the foundation of Einstein's general theory of 
relativity; subsequently, many experimental efforts focused on testing the equiva- 
lence principle in the search for limits of ge nera l relativity. Thus, the early tests of 
the EP were further improved b y Ro ll et alP^ to one part in 10 11 . Most recently, 
a University of Washington group^^l has improved upon Dicke's verification of the 
EP by several orders of magnitude, reporting Mg/Mi — 1 < 1.4 x 10~ 13 . 

The nature of gravity is fundamental to our understanding of our solar system, 
the galaxy and the structure and evolution of the universe. This importance moti- 
vates various precision tests of gravity both in laboratories and in space. To date, 
the experimental evidence for gravitational physics is in agreement with the gen- 
eral theory of relativity; however, there are a number of reasons to question the 
validity of this theory. Despite the success of modern gauge field theories in describ- 
ing the electromagnetic, weak, and strong interactions, it is still not understood 
how gravity should be described at the quantum level. In theories that attempt 
to include gravity, new long-range forces can arise in addition to the Newtonian 
inverse-square law. Even at the purely classical level, and assuming the validity of 
the equivalence principle, Einstein's theory does not provide the most general way 
to establish the space-time metric. Regardless of whether the cosmological constant 
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should be included, there are also important reasons to consider additional fields, 
especially scalar fields. 

Although scalar fields naturally appear in the modern theories, their inclusion 
predicts a non-Einsteinian behavior of gravitating systems. These deviations from 
general relativity lead to a violation of the EP, modification of large-scale grav- 
itational phenomena, and cast doubt upon the constancy of the "constants." In 
particular, the recent work in sca lar-tensor extensions of gravity that are consistent 
with present cosmolog 
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predicts a violation of the EP 
at levels of 10 -13 to 10 -18 . This prediction motivates new searches for very small 
deviations of relativistic gravity from general relativity and provides a robust the- 
oretical paradigm and constructive guidance for further gravity experiments. As a 
result, this theoretical progress has given a new strong motivation for high precision 
tests of relativistic gravity and especially those searching for a possible violation of 
the equivalence principle. Moreover, because of the ever increasing practical sig- 
nificance of the general theory of relativity (i.e. its use in spacecraft navigation, 
time transfer, clock synchronization, standards of time, weight and length, etc) this 
fundamental theory must be tested to increasing accuracy. 

Today Lunar Laser Ranging (LLR) is well positioned to address the challenges 
presented above. The installation of the cornercube retroreflectors on the lunar 
surface more than 35 years ago with the Apollo 11 lunar landing, initiated a unique 
program of lunar laser ranging tests of the EP. LLR provides a set of highly accurate 
distance measurements between an observatory on the Earth and a corner cube 
retroreflector on the Moon which is then used to determine whether, in accordance 
with the EP, these astronomical bodies are both falling towards the Sun at the same 
rate, despite their different masses and compositions. These tests of the EP with 
LLR were among the science goals of the Apollo project. Today this continuing 
legacy of the Apollo program^ constitutes the longest running experiment from 
the Apollo era; it is also the longest on- going experiment in gravitational physics. 

Analyses of laser ranges to the Moon have provided increasingly stringent limits 
on any violation of the EP; they also enabled accurate determinations of a number 
of relativistic gravity parameters. Ranges started in 1969 and have continued with 
a sequence of improvements for 35 years. Data of the last decade are fit with an rms 
residual of 2 cm. This accuracy permits an EP test for the difference in the ratio 
of the gravitational and inertial masses for the Earth and Moon with uncertainty 
of 1.4 x 10~ 13 (see Refs. l9T | II 1 1 1) . The precise LLR data contribute to many areas 
of fundamental and gravitational physics, lunar science, astronomy, and geophysics. 
With a new LLR station in progress and the possibility of new retro-reflectors on 
the Moon, lunar laser ranging remains on the front of gravitational physics research 
in the 21st century. 

This paper focuses on the tests of the EP with LLR. To that extent, Section [2] 
discusses the LLR history, experimental technique, and the current state of the 
effort. Section [3] is devoted to the discussion of the tests of the EP with the Moon. 
It also introduces various "flavors" of the EP and emphasizes the importance of 
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the Earth and Moon as two test bodies to explore the Strong Equivalence Principle 
(SEP). Section |4] describes the existing LLR data including the statistics for the 
stations and reflectors, observational selection effects, and distributions. Section [5] 
introduces and characterizes the modeling and analysis techniques, focusing on the 
tests of the EP. In Section [6] we discuss the details of the scientific data analysis 
using the LLR data set for tests of the EP. We present solutions for the EP and also 
examine the residuals in a search for any systematic signatures. Section[7]focuses on 
the effects derived from the precision tests of the EP. Section [8] introduces the near 
term emerging opportunities and addresses their critical role for the future progress 
in the tests of the equivalence principle with lunar laser ranging. We conclude with 
a summary and outlook. 

2. Lunar Laser Ranging: History and Techniques 

LLR accurately measures the round-trip time of flight for a laser pulse fired from an 
observatory on the Earth, bounced off of a corner cube retroreflector on the Moon, 
and returned to the observatory. The currently available set of LLR measurements 
is more than 35 years long and it has become a major tool to conduct precision tests 
of the EP in the solar system. Notably, if the EP were to be violated this would 
result in an inequality of gravitational and inertial masses and thus, it would lead to 
the Earth and the Moon falling towards the Sun at slightly different rates, thereby 
distorting the lunar orbit. Thus, using the Earth and Moon as astronomical test 
bodies, the LLR experiment searches for an EP-violation-induced perturbation of 
the lunar orbit which could be detected with the available ranges. 

In this Section we discuss the history and current state for this unique experi- 
mental technique used to investigate relativistic gravity in the solar system. 

2.1. Lunar Laser Ranging History 

The idea of using the orbit of the Moon to test foundations of general relativ- 
ity belongs to R. H. Dicke, who in early 1950s suggested using powerful, pulsed 
searchligh ts on the Earth to illuminate corner retroreflectors on the Moon or a 
spacecrafts^ The initial proposal was similar to what today is known as astro- 
metric optical navigation which establishes an accurate trajectory of a spacecraft by 
photographing its position against the stellar background. The progress in quantum 
optics that resulted in the invention of the laser introduced the possibility of ranging 
in early 1960s. Lasers — with their spatial coherence, narrow spectral emission, small 
beam divergence, high power, and well-defined spatial modes — are highly useful for 
many space applications. Precision laser ranging is an excellent example of such a 
practical use. The technique of laser Q-switching enabled laser pulses of only a few 
nanoseconds in length, which allowed highly accurate optical laser ranging. 

Initially the methods of laser ranging to the Moon were analogous to radar 
ranging, with laser pulses bounced off of the lunar surface. A number of these 
early lunar laser ranging experiments were performed in the early 1960's, both 
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at the Massachusetts Institute of Tec hnolo gy and in the former Soviet Union at 
the Crimean astrophysics observatory !^?! However, these lunar surface ranging 
experiments were significantly affected by the rough lunar topography illuminated 
by the laser beam. To overcome this difficulty, deployment of a compact corner 
retroreflector package on the lunar surface was proposed as a part of the unmanned, 
soft-landing Surveyor missions, a proposal that was never realized.^ It was in the 
late 1960's, with the beginning of the NASA Apollo missions, that the concept of 
laser ranging to a lunar corner-cube retroreflector array became a reality. 

The scientific potential of lunar laser ranging led to the placement of retrore- 
flector arrays on the lunar surface by the Apollo astronauts and the unmanned 
Soviet Luna missions to the Moon. The first deployment of such a package on the 
lunar surface took place during the Apollo 1 1 mission (Figure [1]) in the summer of 
1969 and LLR became a reality^! Additional retroreflector packages were set up on 
the lunar surface by the Apollo 14 and 15 astronauts (Figure [2]). Two French-built 
retroreflector arrays were on the Lunokhod 1 and 2 rovers placed on the Moon by 
the Soviet Luna 17 and Luna 21 missions, respectively (Figure [3^). Figure [3}d shows 
the LLR reflector sites on the Moon. 

The first successful lunar laser ranges to the Apollo 11 retroreflector were made 
with the 3.1m telescope at Lick Observatory in northern Californi£0S2lThe ranging 
system at Lick was designed solely for quick acquisition and confirmation, rather 
than for an extended program. Ranges started at the McDonald Observatory in 
1969 shortly after the Apollo 11 mission, while in the Soviet Union a seq uence of 
laser ranges was made from the Crimean astrophysical observatory! ^^ A lunar 
laser ranging program has been carried out in Australia at the Orroral Observa- 
torjQ. Other lunar laser range detections were reported by the Air Force Cambridge 
Research Laboratories Lunar Ranging Observatory in ArizonaP, the Pic du Midi 
Observatory in FrancJl^, and the Tokyo Astronomical ObservatorjI^Ql. 

While some early efforts were brief and demonstrated capability, most of the 
scientific results came from long observing campaigns at several observatories. The 
LLR effort at McDonald Observatory in Texas has been carried out from 1969 to 
the present. The first sequence of observations was made from the 2.7 m telescope. 
In 1985 ranging operations were moved to the McDonald Laser Ranging System 
(MLRS) and in 1988 the MLRS was moved to its present sit<B The MLRS has 
the advantage of a shorter laser pulse and improved range accuracy over the earlier 
2.7 m system, but the pulse energy and aperture are s maller . From 1978 to 1980 a 
set of observations was made from Orroral in Austr alia.™ Accurate observations 
began at the Observatoire de la Cote dAzur (OCA) in and continue to the 

present, though first detections were demonstrated earlier. Ranges were made from 



a The Lick Observatory website: http://www.ucolick.org/ 

b The Orroral Observatory website: http://www.ga.gov.au/nmd/geodesy/slr/index.htm 
c The McDonald Observatory website: http://www.csr.utexas.edu/mlrs/ 
d The Observatoire de la Cote dAzur website: http://www.obs-nice.fr/ 
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Fig. 2. Apollo 14 (left) and Apollo 15 (right) LLR retrorefiector arrays. 



the Haleakala Observatory on the island of Maui in the Hawaiian chain from 1984 
to 199(3. 

Two modern stations which have demonstrated lunar capability are the Wettzell 
Laser Ranging System in German^ and the Matera Laser Ranging Station in 
Italj@. Neither is operational for LLR at present. The Apache Point Obser- 
vatory Lunar Laser ra nging Operation (APOLLO) was recently built in New 
Mexico 62 110 63 64 94 



e The Haleakala Observatory website: http://koa.ifa.hawaii.edu/Lure/ 

f The Wettzell Observatory website: http: //www. wettzell . if ag . de/| 

g The Matera Observatory: http://www.asi.it/html/eng/asicgs/geodynamics/mlro.html 



2, 2009 7:33 WSPC/INSTRUCTION FILE LLR-EP-tests 




The two stations that have produced LLR observations routinely for decades 
are t he Mc Donald Laser Ranging System (MLRSpSI in the United States and the 
station in France. 

2.2. LLR and Fundamental Physics Today 

The analyses of LLR measurements contribute to a wide range of scientific disci- 
plines, and are solely responsible for production of the lunar ephemeris. For a general 
review of LLR see Ref. [39] An independent analysis for Ref. [20] gives geodetic and 
astronomical results. The interior, tidal response , and phy sical librations (rotational 
variations) of the Moon are all probed by LLRj ' making it a valuable tool 
for lunar science. 

The geometry of the Earth, Moon, and orbit is shown in Figure [4] The mean 
distance of the Moon is 385,000 km, but there is considerable variation owing to the 
orbital eccentricity and perturbations due to Sun, planets, and the Earth's J2 zonal 
harmonic. The solar perturbations are thousands of kilometers in size and the lunar 
orbit departs significantly from an ellipse. The sensitivity to the EP comes from the 
accurate knowledge of the lunar orbit. The equatorial radii of the Earth and Moon 
are 6378 km and 1738 km, respectively, so that the lengths and relative orientations 
of the Earth-Moon vector, the station vector, and the retroreflector vector influence 
the range. Thus, not only is there sensitivity of the range to anything which affects 
the orbit, there is also sensitivity to effects at the Earth and Moon. These various 
sensitivities allow the ranges to be analyzed to determine many scientific parameters. 

Concerning fundamental physics, LLR currently provides the most viable solar 
system technique for testing the Strong Equivalence Principle (SEP)-the statement 
that all forms of mass and energy contribute equivalent quantities of inertial and 
gravitational mass (see discussion in the following Section) . The SEP is more restric- 
tive than the weak EP, which applies to non-gravitational mass-energy, effectively 
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probing the compositional dependence of gravitational acceleration. 




Fig. 4. Lunar laser ranging accurately measures the distance between an observatory on Earth 
and a retroreflector on the Moon. 



In addition to the SEP, LLR is capable of measuring the time variation of New- 
ton's gravitational constant, G, providing the strongest limit available for the vari- 
ability of this "constant." LLR can also precisely measure the de Sitter precession- 
effectively a spin-orbit coupling affecting the lunar orbit in the frame co-moving with 
the Earth-Moon system's motion around the Sun. The LLR results are also consis- 
tent with the existence of gravitomagnetism within 0.1% of the predicted levePES 
the lunar orbit is a unique laboratory for gravitational physics where each term in 
the parametrized post-Newtonian (PPN) relativistic equations of motion is verified 
to a very high accuracy. 

A comprehensive paper on tests of gravitational physics is Williams et alff^ A 
recent test of the EP is in Ref. [8] and other general relativity tests are in Ref. 11081 
An overview of the LLR gravitational physics tests is given by Nodtvedt!^ Reviews 
of various tests of relativity, including the contribution by LLR, are given in papers 
by Will .HQ^I 103 | Q ur recen ^ paper, Ref. 1 1 101 describes the model improvements 
needed to achieve the mm-level accuracy for LLR. The most recent LLR results for 
gravitational physics are given in our recent paper of Ref. 11111 



3. Equivalence Principle and the Moon 

Since Newton, the question about equality of inertial and passive gravitational 
masses arises in almost every theory of gravitation. Thus, almost one hundred 
years ago Einstein postulated that not only mechanical laws of motion, but also 
all non-gravitational laws should behave in freely falling frames as if gravity were 
absent. If local gravitational physics is also independent of the more extended grav- 
itational environment, we have what is known as the strong equivalence principle. 
It is this principle that predicts identical accelerations of compositionally different 
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objects in the same gravitational field, and also allows gravity to be viewed as a 
geometrical property of space-time-leading to the general relativistic interpretation 
of gravitation. 

The Equivalence Principle tests can therefore be viewed in two contexts: tests of 
the foundations of the standard model of gravity (i.e. general theory of relativi ty), or 
as searches for new physics because, as emphasized by Damour and colleagues) I I 
almost all extensions to the standard model of particle physics generically predict 
new forces that would show up as apparent violations of the EP. The SEP became 
a foundation of Einstein's general theory of relativity proposed in 1915. Presently, 
LLR is the most viable solar system technique for accurate tests of the SEP, pro- 
viding stringent limits on any possible violation of general relativity - the modern 
standard theory of gravity. 

Below we shall discuss two different "flavors" of the Principle, the weak and the 
strong forms of the EP that are currently tested in various experiments performed 
with laboratory tests masses and with bodies of astronomical sizes. 

3.1. The Weak Form of the Equivalence Principle 

The weak form of the EP (the WEP) states that the gravitational properties of 
strong and electro-weak interactions obey the EP. In this case the relevant test- 
body differences are their fractional nuclear-binding differences, their neutron-to- 
proton ratios, their atomic charges, etc. Furthermore, the equality of gravitational 
and inertial masses implies that different neutral massive test bodies will have the 
same free fall acceleration in an external gravitational field, and therefore in freely 
falling inertial frames the external gravitational field appears only in the form of a 
tidal interactional Apart from these tidal corrections, freely falling bodies behave 
as if external gravity were absent.^ General relativity and other metric theories 
of gravity assume that the WEP is exact. However, extensions of the standard 
model of particle physics that contain new macroscopic-range quantum fields predict 
quantum exchange forces that generically violate the WEP becau se the y couple to 
generalized "charges" rather than to mass/energy as does gravityPZEEl 

In a laboratory, precise tests of the EP can be made by comparing the free fall 
accelerations, a\ and 02, of different test bodies. When the bodies are at the same 
distance from the source of the gravity, the expression for the equivalence principle 
takes an elegant form: 



where Mg and Mj represent gravitational and inertial masses of each body. The 
sensitivity of the EP test is determined by the precision of the differential accel- 
eration measurement divided by the degree to which the test bodies differ (e.g. 
composition). 

Since the early days of general relativity, Einstein's version of the Equivalence 
Principle became a primary focus of many experimental efforts. Various experiments 




(1) 
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have been performed to measure the ratios of gravitational to inertial masses of 
bodies. Recent experiments on bodies of laboratory dimensions verify the WEP 
to a fractional precision A(M G /M>) < 10~ n by Roll et alP^ to < 10" 12 by 
Refs. [T3l HI] and more recently to a precision of < 1.4 x 10~ 13 in Ref. [5] The 
accuracy of these experiments is sufficiently high to confirm that the strong, weak, 
and electromagnetic interactions each contribute equally to the passive gravitational 
and inertial masses of the laboratory bodies. 

This impressive evidence for laboratory bodies is incomplete for astronomical 
body scales. The experiments searching for WEP violations are conducted in lab- 
oratory environments that utilize test masses with negligible amounts of gravita- 
tional self-energy and therefore a large scale experiment is needed to test the postu- 
lated equality of gravitational self-energy contributions to the inertial and passive 
gravitational masses of the bodied^. Once the self-gravity of the test bodies is 
non-negligible (currently with bodies of astronomical sizes only), the corresponding 
experiment will be testing the ultimate version of the EP - the strong equivalence 
principle, that is discussed below. 



3.2. The Strong Form of the Equivalence Principle 

In its strong form the EP is extended to cover the gravitational properties result- 
ing from gravitational energy itself. In other words, it is an assumption about the 
way that gravity begets gravity, i.e. about the non-linear property of gravitation. 
Although general relativity assumes that the SEP is exact, alternate metric theo- 
ries of gravity such as those invol ving scalar fi elds, and other extensions of gravity 
theory, typically violate the SEP.ESmSEZES! p or the SEP case, the relevant test 
body differences are the fractional contributions to their masses by gravitational 
self-energy. Because of the extreme weakness of gravity, SEP test bodies that differ 
significantly must have astronomical sizes. Currently, the Earth-Moon-Sun system 
provides the best solar system arena for testing the SEP. 

A wide class of metric theories of gravity are described by the parametrized 
post-Newtonian which allows one to describe within a common 

framework the motion of celestial bodies in external gravitational fields. Over the 
last 35 years, the PPN formalism has become a useful framework for testing the SEP 
for extended bodies. To facilitate investigation of a possible violation of the SEP, 
in that f ormali sm the ratio between gravitational and inertial masses, Mq/Mj, is 
expresse d^lMi as 



Mi 



SEP 

where M is the mass of a body, U is the body's gravitational self-energy (U < 

0), Mc 2 is its to tal mass-energy, and r\ is a dimensionless constant for SEP 
violation!§SIIMII 

Any SEP violation is quantifie d by the parameter r\. In fully-conservative, 
Lorentz-invariant theories of gravity ^ ^ I IO^ I the SEP parameter is related to the 
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(3) 




(4) 



For a sphere with a radius R and uniform density, U/Mc 2 — —ZGM/hRc 2 = 
—3v E /10c 2 , where ve is the escape velocity. Accurate evaluation for solar system 
bodies requires numerical integration of the expression of Eq. ((4]). Evaluating the 
standard solar models results in (U/Mc 2 )s ~ —3.52 x 10~ 6 . Because gravitational 
self-energy is proportional to M 2 (i.e. U/Mc 2 ~ M) and also because of the extreme 
weakness of gravity, the typical values for the ratio (U/Mc 2 ) are ~ 10 -25 for bod- 
ies of laboratory sizes. Therefore, the experimental accuracy of a part in 10 13 (see 
Ref. [5]) which is so useful for the WEP is not a useful test of how gravitational self- 
energy contributes to the inertial and gravitational masses of small bodies. To test 
the SEP one must utilize planetary-sized extended bodies where the ratio Eq. (j4|) 
is considerab ly higher. 

Nordtvedl ! 1 1 1 suggested several solar system experiments for testing the 
SEP. One of these was the lunar test. Another, a search for the SEP effect in the 
motion of the Trojan asteroids, was carried out by Orellana and Vucetich F^ I T^ I 
Interplanetary spacecraft tests have been considered by Anderson et alP and dis- 
cussed by Anderson and Williams.^ An experime nt e mploying existing binary pulsar 
data has been proposed by Damour and SchaferPH It was pointed out that binary 
pulsars may provide an ex cellent possibility for testing the SEP in the new regime 
of strong self-gr ayity j^ ^ l however the corresponding tests have yet to reach com- 
petitive accuracySSIMl To date, the Earth-Moon-Sun system has provided the most 
accurate test of the SEP with LLR being the available technique. 

3.3. Equivalence Principle and the Earth-Moon system 

The Earth and Moon are large enough to have significant gravitational self energies 
and a lunar test of the equivalence principle was proposed by Nordtvedt.^ Both 
bodies have differences in their compositions and self energie s and the Sun provides 
the external gravitational acceleration. For the EarthSSHnSl a nume rical evaluation 
of Eq. (0]) yields: 



The two evaluations, with different Earth models, differ by only 0.1%. (A uniform 
Earth approximation is 10% smaller in magnitude.) A Moon model, with an iron 
core ~20% of its radius, gives 




(5) 




(6) 
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The subscripts E and M denote the Earth and Moon, respectively. The lunar value 
is only 1% different from the uniform density approximation which demonstrates its 
insensitivity to the model. The lunar value was truncated to two digits in Ref. 11051 
For the SEP effect on the Moon's position with respect to the Earth it is the 
difference of the two accelerations and self-energy values which is of interest. 



The Jet Propulsion Laboratory's (JPL) program which integrates the orbits 
of the Moon and planets considers accelerations due to Newtonian, geophysical 
and post-Newtonian effects. Considering just the modification of the point mass 
Newtonian terms, the equivalence principle enters the acceleration aj of body j as 



Mc 2 J 



G TFT, >>^- W 



j k 3 k 



where Tju = Yk — tj is the vector from accelerated body j to attracting body k and 
Tjk — \rjk\. For a more through discussion of the integration model see Ref. [571 

The dynamics of the three-body Sun-Earth-Moon system in the solar system 
barycentric inertial frame provides the main LLR sensitivity for a possible violation 
of the equivalence principle. In this frame, the quasi-Newtonian acceleration of the 
Moon with respect to the Earth, a = aj,; — sle, is calculated to be: 

* r EM (M G \ Y S M . ( M G \ Y SE 



where [i* = He(Mg/Mi)m + Hm{Mg/Mi)e and fj,k = GMk- The first term on the 
right-hand side of Eq.©, is the acceleration between the Earth and Moon with the 
remaining pair being the tidal acceleration expression due to the solar gravity. The 
above acceleration is useful for either the weak or strong forms of the EP. 

For the SEP case, rj enters when expression Eq. ((2|) is is combined with Eq. , 



a = -/j, — h^s 

r EM 



r, 



\Mc 2 ) E r% E \Mc 2 ) M r% M 



(10) 



SE ' SM . 

The presence of r] in //* modifies Kepler's third law to n 2 a 3 = [i* for the relation 
between semimajor axis a and mean motion n in the elliptical orbit approximation. 
This term is notable, but in the LLR solutions /j.e + I^m is a solution parameter, or 
at least uncertain (see Sec. [7J, so this term does not provide a sensitive test of the 
equivalence principle, though its effect is implicit in the LLR solutions. The second 
term on the right-hand side with the differential acceleration toward the Sun is the 
Newtonian tidal acceleration. The third term involving the self energies gives the 
main sensitivity of the LLR test of the equivalence principle. Since the distance to 
the Sun is ^390 times the distance between the Earth and Moon, the last term, is 
approximately rj times the difference in the self energies of the two bodies times the 
Sun's acceleration of the Earth-Moon center of mass. 
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Treating the EP related tidal term as a perturbation NordtvedlP^ found a po- 
larization of the Moon's orbit in the direction of the Sun with a radial perturbation 



Ar = S 



(Mg\ _ (Ma 
\Mj) E \Mj 



cosD, (11) 



MJ 

where S is a scaling factor of about —2.9 x 10 13 mm (see Refs. \71 \ l32 | l33|) . For the 
SEP, combining Eqs. (J2]) and yields 

U E U m 



Ar = Srj 



cosD, (12) 



M E c 2 M M c 2 
Ar = C r) cos D. (13) 

Applying the difference in numerical values for self-energy for the Earth and Moon 
Eq. gives a value of Co of about 13 m (see Refs. 175 } I32ll33|) . In general relativity 
r\ = 0. A unit value for r\ would produce a displacement of the lunar orbit about 
the Earth, causing a 13 m monthly range modulation. See subsection 16.11 for a 
comparison of the theoretical values of S and Co with numerical results. This effect 
can be generalized to all similar three body situations. 

In essence, LLR tests of the EP compare the free-fall accelerations of the Earth 
and Moon toward the Sun. Lunar laser-ranging measures the time-of-flight of a laser 
pulse fired from an observatory on the Earth, bounced off of a retroreflector on the 
Moon, and returned to the observatory (see Refs. I3"9l ITTj) . If the Equivalence Princi- 
ple is violated, the lunar orbit will be displaced along the Earth-Sun line, producing 
a range signature having a 29.53 day synodic period (different from the lunar orbit 
period of 27 days). The first LLR tests of the EP wer e published in 1976 (see 11041 
ESI). Since then the precision of the test has increasecP E 39 ! ! 7 ! 105 ! 106 ! 61 ! 8 ! 108 ! !!^ 
until modern results are improved by two orders-of-magnitude. 

3.4. Equivalence Principle and Acceleration by Dark Matter 

At the scales of galaxies and larger there is evidence for unseen dark matter. Thus, 
observations of disk galaxies imply that the circular speeds are approximately in- 
dependent of distance to the center of the galaxy at large distances. The standard 
explanation is that this is due to halos of unseen matter that makes up around 90% 
of the total mass of the galaxies!^ The same pattern repeats itself on larger and 
larger scales, until we reach the cosmic scales where a baryonic density compatible 
with successful big bang nucleosynthesis is less tha n 10% of the density predicted 
by inflation, i.e. the critical density. Braginsky et alP^ES have studied the effect of 
dark matter bound in the galaxy but unbound to the solar system. Such galactic 
dark matter would produce an anisotropy in the gravitational background of the 
solar system. 

A possible influence of dark matter on the Earth-Moon system has been consid- 
ered by NordtvedtP^, who has pointed out that LLR can also test ordinary matter 
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interacting with galactic dark matter. It was suggested that LLR data can be used 
to set experimental limits to the density of dark matter in the solar system by 
studying its effect upon the motion of the Earth-Moon system. The period of the 
range signature is the sidereal month, 27.32 days. An anomalous acceleration of 
10 -15 m/s 2 would cause a 2.5 cm range perturbation. At this period there are also 
signatures due to other solution parameters: one component of station location, 
obliquity, and orbital mean longitude. These parameters are separable because they 
contribute at other periods as well, but they are complications to the dark matter 
test. 

In 1995, Nordtvedt Miiller, and Soffel published an upper limit of 3x 10~ 16 m/s 2 
for a possible differential acceleration in coupling of dark matter to the different 
compositions of Earth and Moon. This represented a stronger constraint by a factor 
of 150 than was achieved by the laboratory experiments searching for differential 
cosmic ac celeration rates between beryllium and copper and between beryllium and 

4. Data 

The accuracy and span of the ranges limit the accuracy of fit parameters. This 
section describes the data set that is used to perform tests of the Equivalence 
Principle with LLR. The data taking is a day-to-day operation at the McDonald 
Laser Ranging System (MLRS) and the Observatoire de la Cote d'Azur (OCA) 
stations. 

LLR has remained a viable experiment with fresh results over 35 years because 
the data accuracies have improved by an order of magnitude. See Section 4.1 below 
for a discussion and illustration (Figure [5]) of that improvement. The International 
Laser Ranging Service (ILRsjEI provides lunar laser ranging data and their related 
products to support geodetic and geophysical research activities. 

4.1. Station and Reflector Statistics 

LLR data have been acquired from 1969 to the present. Each measurement is the 
round-trip travel time of a laser pulse between a terrestrial observatory and one of 
four corner cube retroreflectors on the lunar surface. A normal point is the standard 
form of an LLR datum used in the analysis. It is the result of a statistical combining 
of the observed transit times of several individual photons detected by the observing 
instrument within a relatively short time, typically a few minutes to a few tens of 
minutes. 

T he currently operating LLR stations , M cDonald Laser Ranging System in 
TexaJ^and Observatoire de la Cote d'AzuJ^, typically detect 0.01 return photons 
per pulse during normal operation. A typical "normal point " is constructed from 
3-100 return photons, spanning 10-45 minutes of observation!^ 



h International Laser Ranging Service (ILRS) website at http://ilrs.gsfc.nasa.gov/index.html 
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Historical Accuracy of LLR Data 



35 




I 1 1 1 



1970 1975 1980 1985 1990 1995 2000 2005 

Fig. 5. Annual rms residuals of LLR data from 1970 to 2004. 

The LLR data set for analysis has observations from McDonald Observatory, 
Haleakala Observatory, and OCA. Figure [5] shows the weighted RMS residual for 
each year. Early accuracies using the McDonald Observatory's 2.7 m telescope hov- 
ered around 25 cm. Equipment improvements decreased the ranging uncertainty to 
~15 cm later in the 1970s. In 1985 the 2.7 m ranging system was replaced with the 
MLRS. In the 1980s lunar ranges were also received from Haleakala Observatory on 
the island of Maui, Hawaii, and OCA in France. Haleakala ceased lunar operations 
in 1990. A sequence of technical improvements decreased the rms residual to the 
current ~2 cm of the past decade. The 2.7 m telescope had a greater light gathering 
capability than the newer smaller aperture systems, but the newer systems fired 
more frequently and had a much improved range accuracy. The new systems do not 
distinguish returning photons against the bright background near full Moon, which 
the 2.7 m telescope could do. There are some modern eclipse observations. 

The first LLR test of the EP used 1523 normal points up to May 1975 with 
accuracies of 25 cm. By April 2004, the data set has now grown to more than 
15,554 normal points spanning 35 years, and the recent data is fit with ~2 cm rms 
scatter. Over time the post-fit rms residual has decreased due to improvements at 
both the McDonald and the OCA sites. Averaged over the past four years there 
have been a total of several hundred normal points per year. 

The full LLR data set is dominated by three stations: the McDonald Station 
in Texas, the OCA station at Grasse, France, and the Haleakala station on Maui, 
Hawaii. At present, routine ranges are being obtained only by the MLRS and OCA. 
Figure[3j3 shows the distribution of the lunar retroreflectors. Over the full data span 
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Fig. 6. Corner-cube prisms are optical devices that return any incident light back in exactly the 
direction from which it came (left). An array of corner-cubes makes up the Apollo 11 lunar laser 
reflector (right). 



78% of the ranges come from Apollo 15, 10% from Apollo 11, 9% from Apollo 14, 
3% from Lunokhod 2, and nothing from Lunokhod 1 (lost). 

The notable improvement of the LLR data set with time implies comparable 
improvement in the determination of the solution parameters. Data from multiple 
ranging sites to multiple retrorefiectors are needed for a robust analysis effort. 

4.2. Observational Influences and Selection Effects 

To range the Moon the observatories on the Earth fire a short laser pulse toward the 
target retroreflector array. The outgoing laser beam is narrow and the illuminated 
spot on the Moon is a few kilometers across. The retrorefiectors are made up of 
arrays of corner cubes: 100 for Apollos 11 and 14, 300 for Apollo 15, and 14 for 
the Lunokhods. At each corner cube (Figure [6]) the laser beam enters the front face 
and bounces off of each of the three orthogonal faces at the rear of the corner cube. 
The triply reflected pulse exits the front face and returns in a direction opposite 
to its approach. The returning pulse illuminates an area around the observatory 
which is a few tens of kilometers in diameter. The observatory has a very sensitive 
detector which records single photon arrivals. Color and spatial filters are used to 
eliminate much of the background light. Photons from different laser pulses have 
similar residuals with respect to the expected round-trip time of flight and are 
thus separated from the widely scattered randomly arriving background photons. 
The resulting "range" normal point is the round trip light time for a particular 
firing time. (For more details on satellite and lunar laser ran ging instr umentation, 
experimental set-up. and operations, consult papers by Degnan^ 35 36 anf | g amam 
et alPl 

The signal returning from the Moon is so weak that single photons must be de- 
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tected. Not all ranging attempts are successful and the likelihood of success depends 
on the conditions of observation. Observational effects may influence the strength 
of the signal, the background light which competes with the detection of the re- 
turning laser signal, the width of the outgoing or returning beam, and the telescope 
pointing. Some of these observational influences select randomly and some select 
systematically, e.g. with phase of Moon, time of day, or time of year. Selection with 
phase influences the equivalence principle test. This subsection briefly discusses 
these observational influences and selection effects. 

The narrow laser beam must be accurately pointed at the target. Seeing, a mea- 
sure of the chaotic blurring of a point source during the transmission of light through 
the atmosphere, affects both the outgoing laser beam and the returning signal. The 
beam's angular spread, typically a few seconds of arc ("), depends on atmospheric 
seeing so the spot size on the Moon is a few kilometers across at the Moon's distance 
(use 1.87 km/"). The amount of energy falling on the retroreflector array depends 
inversely on that spot area. At the telescope's detector both a diaphragm restricting 
the field of view and a (few Angstrom) narrow-band color filter reduce background 
light. When the background light is high the diaphragm should be small to reduce 
the interference and increase the signal-to-noise ratio. When the seeing is poor the 
image size increases and this requires a larger diaphragm. 

The phase of the Moon determines whether a target retroreflector array is il- 
luminated by sunlight or is in the dark. These phase effects include the following 
influences. 

a) The target illumination determines the amount of sunlight scattered back 
toward the observatory from the lunar surface near the target. A sunlit sur- 
face increases the noise photons at the observatory's detector and decreases 
the signal to noise ratio. 

b) The pointing technique depends on solar illumination around the target 
array. Visual pointing is used when the target is sunlit while more difficult 
offset pointing, alignment using a displaced illuminated feature, is used 
when the target is dark. 

c) Retroreflector illumination by sunlight determines solar heating of the ar- 
ray and thermal effects on the retroreflector corner cubes. A thermal gra- 
dient across a corner cube distorts the optical quality and spreads the re- 
turn beam. The Lunokhod corner cubes are about twice the size of the 
Apollo corner cubes and are thus more sensitive to thermal effects. Also, 
the Lunokhod corner cubes have a reflecting coating on the three reflecting 
back sides while the Apollo corner cubes depend on total internal reflection. 
The coating improves the reflected strength for beams that enter the front 
surface at an angle to the normal, where the Apollo efficiency decreases, 
but it also heats when sunlit. Thus, the Lunokhod arrays have greater ther- 
mal sensitivity and are more difficult targets when heated by sunlight. A 
retroreflector in the dark is in a favorable thermal environment, but the 
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telescope pointing is more difficult. 

Whether the observatory is experiencing daylight or night determines whether 
sunlight is scattered toward the detector by the atmosphere. As the Moon's phase 
approaches new, the fraction of time the Moon spends in the observatory's daylight 
sky increases while the maximum elevation of the Moon in the night sky decreases, 
so atmosphere-scattered sunlight is correlated with lunar phase. 

The beam returning from the Moon cannot be narrower than the diffraction 
pattern for a corner cube. The diffraction pattern of a corner cube has a six-fold 
shape that depends on the six combinations of ways that light can bounce off of 
the three orthogonal reflecting faces. An approximate computation for green laser 
light (0.53 jLtm) gives 7 arcsec for the angular diameter of an Airy diffraction disk. 
The larger Lunokhod corner cubes would give half that diffraction pattern size. 
Thermal distortions, imperfections, and contaminating dust can make the size of 
the returning beam larger than the diffraction pattern. So the returning spot size 
on the Earth is ~30 km across for green laser light. The power received by the 
telescope depends directly on the telescope's collecting area and inversely on the 
returning spot area. Velocity-caused aberration of the returning beam is roughly 1" 
and is not a limitation since it is much smaller than the diffraction pattern. 

There are geometrical selection effects. For the two operational northern rang- 
ing stations the Moon spends more time above the horizon when it is at northern 
declinations and less when south. Also, atmospheric effects such as seeing and ab- 
sorption increase at low elevation. Consequently, there is selection by declination of 
the Moon. This, along with climate, causes seasonal selection effects. A station can 
only range the Moon when it is above the horizon which imposes selection at the 
24 hr 50.47 min mean interval between meridian crossings. 

The best conditions for ranging occur with the Moon located high in a dark sky 
on a night of good seeing. A daylight sky adds difficulty and full Moon is even more 
difficult. A retroreflector in the dark benefits from not being heated by the Sun, 
but aiming the laser beam is more difficult. New Moon occurs in the daylight sky 
near the Sun and ranging is not attempted since sensitive detectors are vulnerable 
to damage from bright light. 

4.3. Data Distributions 

Observational selection effects shape the data distribution. Several selection effects 
depend on the phase of the Moon and there is a dramatic influence on the distri- 
bution of the number of observations with phase. The elongation of the Moon from 
the Sun is approximated with the angle D, the smooth polynomial representation 
of the difference in the mean longitudes for Sun and Moon. Zero is near new Moon, 
90° is near first quarter, 180° is near full Moon, and 270° is near last quarter. Fig- 
ure^ illustrates the distribution of observations for the decade from 1995-2004 with 
respect to the angle D. The shape of the curve results from the various selection 
effects discussed above. There are no ranges near new Moon and few ranges near full 
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Moon. The currently operating observatories only attempt full Moon ranges during 
eclipses. The original 2.7 m McDonald ranging system transmitted more energy in 
its longer pulse than currently operating systems, which gave it a higher single shot 
signal to noise ratio against a bright background. It could range during full Moon 
as the distribution of the full data set for 1970-2004 shows in Figure [8^,. 

Factors such as weather and the northern hemisphere location of the operating 
stations cause seasonal selection effects. The distribution of the number of observa- 
tions vs the mean anomaly of the Earth-Moon system about the Sun is shown in 
Figure [5^. The annual mean anomaly is zero in the first week of January so that 
the mean anomaly is offset from calendar day of the year by only a few days. There 
is considerable variation in the frequency of observation; the distribution is at its 
highest in fall and winter and at its lowest in summer. 

Other selection effects such as distance and declination also influence the data 
distribution and can be seen with appropriate histograms. Nonuniform data distri- 
butions are one contribution to correlations between solution parameters. 

5. Modeling 

Lunar Laser Ranging measures the range from an observatory on the Earth to a 
retroreflector on the Moon. The center-to-center distance of the Moon from the 
Earth, with mean value 385,000 km, is variable due to such things as orbit eccen- 
tricity, the attraction of the Sun, planets, and the Earth's bulge, and relativistic 
corrections. In addition to the lunar orbit, the range from an observatory on the 
Earth to a retroreflector on the Moon depends on the positions in space of the rang- 
ing observatory and the targeted lunar retroreflector. Thus, the orientation of the 
rotation axes and the rotation angles for both bodies are important. Tidal distor- 
tions, plate motion, and relativistic transformations also come into play. To extract 
the scientific information of interest, it is necessary to accurately model a variety 
of effects. 

The sensitivity to the equivalence principle is through the orbital dynamics. The 
successful analysis of LLR data requires attention to geophysical and rotational 
effects for the Earth and the Moon in addition to the orbital effects. Modeling is 
central to the data analysis. The existing model formulation, and its computational 
realization in computer code, is the product of much effort. This section gives an 
overview of the elements included in the present model. 

5.1. Range Model 

The time-of- flight ( "range" ) calculation consists of the round-trip "light time" from 
a ranging site on the Earth to a retroreflector on the Moon and back to the ranging 
site. This time of flight is about 2.5 sec. The vector equation for the one-way range 
vector p is 



P — T — R s tn + R r fl 7 



(14) 
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Fig. 7. (a) Distribution of last decade of data vs argument D, which has a 29.53 day period, (b) 
Weighted average residual vs D for last decade. 

where r is the vector from the center of the Earth to the center of the Moon, R st n 
is the vector from the center of the Earth to the ranging site, and R r fi is the 
vector from the center of the Moon to the retroreflector array (see Figure |4] for 
more details). The total time of flight is the sum of the transmit and receive paths 
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Fig. 8. (a) Distribution of all data vs argument D. (b) Weighted average residual vs D for all 
data. 

plus delays due to atmosphere and relativistic gravitational delay 



360 



*3 - ti = (pi2 + P23,)/c + Ai atm + At grav . (15) 
The times at the Earth are transmit (1) and receive (3), while the bounce time (2) is 
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Fig. 9. (a) Annual mean anomaly distribution for the last decade of data, (b) Weighted average 
residual vs annual mean anomaly. 



at the Moon. Due to the motion of the bodies the light-time computation is iterated 
for both the transmit and receive legs. Since most effects effectively get doubled, it 
is convenient to think of 1 nsec in the round-trip time as being equivalent to 15 cm 
in the one-way distance. 
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The center of mass of the solar system is treated as unaccelerated. This solar 
system barycenter (SSB) is the coordinate frame for evaluating the above equa- 
tions including relativistic computations. First, the transmit tim e at the station is 
transformed to the SSB coordinate time (called T e ph by StandishP^I approximated 
by TDB), the basic computations are made in that SSB frame, and the computed 
receive time is transformed back to the station's time. 



(*3 - *l)sta = *3 - *1 + Attrana- (16) 

The form of Eq. (fT4|) separates the modeling problem into aspects related to the 
orbit, the Earth, and the Moon. Eq. (|15p shows that time delays must be added 
and Eq. (| 16[) demonstrates modification of the round-trip-time-delay due to choice 
of reference frame. For the discussion below we make a similar separation. The 
dynamics of the orbits and lunar rotation come from a numerical integration, and 
those are the first two topics. Earth and Moon related computations are discussed 
next. The last topic is time delays and transformations. 



5.1.1. Orbit Dynamics, r 

The lunar and planetary orbits and the lunar rotation result from a simultaneous 
numerical integration of the differential equations of motion. The numerical inte- 
gration model is detailed by Standish and WilliamJ^. Ephemerides of the Moon 
and planets plus lunar rotation are available at the Jet Propulsion Laboratory web 



site http: //ssd.j pi .nasa.gov/ 



The numerical integration of the motion of the Moon, planets, and Sun generates 
positions and velocities vs time. The existing model for accelerations accounts for: 

• Newtonian and relativistic point mass gravitational interaction between 
the Sun, Moon, and nine planets. Input parameters include masses, orbit 
initial conditions, PPN parameters j3 and 7, G, and equivalence principle 
parameters (Mq/Mj). 

• Newtonian attraction of the largest asteroids. 

• Newtonian attraction between point mass bodies and bodies with gravita- 
tional harmonics: Earth (J2, J3, J4), Moon (second- through fourth-degree 
spherical harmonics), and Sun (J2V 

• Attraction from tides on both Earth and Moon includes both elastic and 
dissipative components. There is a terrestrial Love number &2 and a time 
delay for each of three frequency bands: semidiurnal, diurnal, and long 
period. The Moon has a different Love number ki and time delay. 



5.1.2. Lunar Rotation Dynamics 

The numerical integration of the rotation of the Moon generates three Euler angles 
and three angular velocities. The torque model accounts for: 
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• Torques from the point mass attraction of Earth, Sun, Venus, Mars and 
Jupiter. The lunar gravity field includes second- through fourth-degree 
terms. 

• Figure-figure torques between Earth (J2) and Moon (J2 and 622)- 

• Torques from tides raised on the Moon include elastic and dissipative com- 
ponents. The formulation uses a lunar Love number ki and time delay. 

• The fluid core of the Moon is considered to rotate separately from the 
mantle. A diss ipati ve torque at the lunar solid-mantle/fluid-core interface 
couples the twd ^^ l There is a coupling parameter and the rotations of both 
mantle and core are integrated. 

• An oblate fluid-core/solid-mantle boundary generates a torque from the 
flow of the fluid along the boundary. This is a recent addition. 



5.1.3. Effects at Earth, R stn 

• The ranging station coordinates include rates for horizontal plate motion 
and vertical motion. 

• The solid-body tides are raised by Moon and Sun and tidal displacements 
on the Earth are scaled by terrestrial Love numbers ft 2 and I2 ■ There is also 
a core-flattening correction for a nearly diurnal term and a "pole tide" due 
to the time- varying part of the spin distortion. 

• The orientation of the Earth's rotation axis includes precession and nuta- 
tion. The body polar (z) axis is displaced from the rotation axis by polar 
motion. The daily rotation includes UT1 variations. A rotation matrix be- 
tween the space and body frames incorporates these effects. 

• The motion of the Earth with respect to the solar system barycenter re- 
quires a Lorentz contraction for the position of the geocentric ranging sta- 
tion. 

A compilation of Earth-related effects has been collected by McCarthy and PetitP^. 



5.1.4. Effects at the Moon, R r fi 

• The Moon-centered coordinates of the retroreflectors are adjusted for solid- 
body tidal displacements on the Moon. Tides raised by Earth and Sun are 
scaled by the lunar displacement Love numbers /12 and li- 

• The rotation matrix between the space and lunar body frames depends 
on the three Euler angles that come from the numerical integration of the 
Euler equations. 

• The motion of the Moon with respect to the solar system barycenter re- 
quires a Lorentz contraction for the position of the Moon-centered reflector. 
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5.1.5. Time Delays and Transformations 

• Atmospheric time delay Ai a tm follows Ref. 57. It includes corrections for 
surface pressure, temperature and humidity which are measured at the 
ranging site. 

• The relativistic time transformation has time- varying terms due to the mo- 
tion of the Earth's center with respect to the solar system barycenter. In 
addition, the displacement of the ranging station from the center of the 
Earth contributes to the time transformation. The transformation changes 
during the ~2.5 sec round-trip time and must be computed for both trans- 
mit and receive times. 

• The propagation of light in the gravity fields of the Sun and Earth causes 
a relativistic time delay Ai grav . 



5.1.6. Fit Parameters & Partial Derivatives 

For each solution parameter in the least-squares fit there must be a partial derivative 
of the "range" with respect to that parameter. The partial derivatives may be 
separated into two types - geometrical and dynamical. 

Geometrical partials of range are explicit in the model for the time of flight. 
Examples are partial derivatives of range with respect to geocentric ranging sta- 
tion coordinates, Moon-centered reflector coordinates, station rates due to plate 
motion, tidal displacement Love numbers hi and I2 for Earth and Moon, selected 
nutation coefficients, diurnal and semidiurnal UT1 coefficients, angles and rates for 
the Earth's orientation in space, and ranging biases. 

Dynamical partials of lunar orbit and rotation are with respect to parameters 
that enter into the model for numerical integration of the orbits and lunar rotation. 
Examples are dynamical partial derivatives with respect to the masses and orbit ini- 
tial conditions for the Moon and planets, the masses of several asteroids, the initial 
conditions for the rotation of both the lunar mantle and fluid core, Earth and Moon 
tidal gravity parameters (&2 and time delay) , lunar moment of inertia combinations 
(B — A)/C and (C — A)/B, lunar third-degree gravity field coefficients, a lunar core- 
mantle coupling parameter, equivalence principle Mq/Mi, PPN parameters (3 and 
7, geodetic precession, solar J2, and a rate of change for the gravitational constant 
G. Dynamical partial derivatives for the lunar and planetary orbits and the lunar 
rotation are created by numerical integration. 

Considering Eqs. (TH| and (fT5|) . the partial derivative of the scalar range p with 
respect to some parameter p takes the form 



dp dp 

d-p = ip -^ (17) 

where p = p/p is the unit vector. From the three terms in Eq. (|14[) . dp/dp depends 
on dr/dp, —dR s tn/dp, and dR r fi/dp. 
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The dynamical partial derivatives of the orbit are represented by dr/dp. Rotation 
matrices are used to transform R stn and R r fi between body and space-oriented 
coordinates, so the partial derivatives of the rotation matrices depend on fit pa- 
rameters involving the Earth and Moon Euler angles such as the Earth rotation, 
precession and nutation quantities and numerous lunar parameters which are sensi- 
tive through the rotation. Only geometrical partials contribute to dTl stn /dp. Both 
dynamical and geometrical partials affect 9R rf i/<9p. 

5.1.7. Computation 

The analytical model has its computational realization in a sequence of computer 
programs. Briefly these programs perform the following tasks. 

a) Numerically integrate the lunar and planetary orbits along with lunar ro- 
tation. 

b) Numerically integrate the dynamical partial derivatives. 

c) Compute the model range for each data point, form the pre- fit residual, and 
compute range partial derivatives. At the time of the range calculation a file 
of integrated partial derivatives for orbits and lunar rotation with respect to 
dynamical solution parameters is read and converted to partial derivatives 
for range with respect to those parameters following Eq. (fT8|) . The partial 
derivative for PPN 7 has both dynamical and geometrical components. 

d) Solve the weighted least-squares equations to get new values for the fit 
parameters. 

e) Generate and plot post-fit residuals. 

A variety of solutions can be made using different combinations of fit parameters. 
Linear constraints between solution parameters can also be imposed. The dynamical 
parameters from a solution can be used to start a new integration followed by new 
fits. The highest quality ephemerides are produced by iterating the integration and 
fit cycle. 

5.1.8. Data Weighting 

A range normal point is composed of from 3 to 100 single photon detections. As the 
normal point comes from the station, the uncertainty depends on the calibration 
uncertainty and the time spread of the detected returned pulse. The latter depends 
on the length of the outgoing laser pulse, spread at the lunar retroreflector due to 
tilt of the array, and detector uncertainty. Gathering more photons reduces these 
return pulse length contributions to the normal point. The analyst can also adjust 
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the weightings according to experience with the residuals. The analysis program 
includes uncertainty associated with the input UT1 and polar motion variations. 



5.1.9. Solar Radiation Pressure 

Solar radiation pressure, like the acceleration from an equivalence principle viola- 
tion, is aligned with the direction from the Sun and it produces a perturbation with 
the 29.53 d synodic period. Thus, this force on the Earth and Moon deserves special 
consideration for the most accurate tests of the equivalence principle. This accel- 
eration is not currently modeled in the JPL software. Here we rely on the analysis 
of VokrouhlickjE3 who considered incident and reflected radiation for both bodies 
plus thermal radiation from the Moon. He finds a solar radiation perturbation of 
—3.65 ± 0.08 mm cosD in the radial coordinate. 



5.1.10. Thermal Expansion 

The peak to peak variation of surface temperature at low latitudes on the Moon is 
nearly 300° C. The lunar "day" is 29.53 days long. This is the same period as the 
largest equivalence principle term so a systematic effect from thermal expansion is 
indicated. The phase of the thermal cycle depends on the retroreflector longitude. 

The Apollo retroreflector arrays and the Lunokhod 1 vehicle with the attached 
retroreflector array are shown in Figure[T][3] The Apollo 11, 14, and 15 retroreflector 
arrays are close to the lunar surface and the center of each array front face is about 
0.3, 0.2, and 0.3 m above the surface, respectively. The Apollo corner cubes are 
mounted in an aluminum plate. The thermal expansion coefficient for aluminum is 
about 2 x 10~ 5 /°C. If the Apollo arrays share the same temperature variations as 
the surface, then the total variation of thermal expansion will be 1 to 2 mm. The 
Lunokhod 2 vehicle is 1.35 m high. From images the retroreflector array appears to 
be just below the top and it is located in front of the main body of the Lunokhod. 
We do not know the precise array position or the thermal expansion coefficient of 
the rover, but assuming the latter is in the range of 1 x 10~ 5 /°C to 3 x 10~ 5 /°C 
then the peak vertical thermal variation will be in the range of 3 to 10 mm. The 
horizontal displacement from the center of the Lunokhod is poorly known, but it 
appears to be ~1 m and the horizontal thermal variation will be similar in size to the 
vertical variation. The thermal expansion cycle is not currently modeled. For future 
analyses, it appears to be possible to model the thermal expansion of the Apollo 
arrays without solution parameters, but a solution parameter for the Lunokhod 2 
thermal cycle expansion seems to be indicated. 

The soil is heated and subject to thermal expansion, but it is very insulating and 
the "daily" thermal variation decreases rapidly with depth. So less displacement is 
expected from the thermal expansion of the soil than from the retroreflector array. 
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6. Data Analysis 

This section presents analysis of the lunar laser ranging data to test the equivalence 
principle. To check consistency, more than one solution is presented. Solutions are 
made with two different equivalence principle parameters and different ways of 
establishing the masses of Earth and Moon. Also, spectra of the residuals after fits, 
the post-fit residuals, are examined for systematics. 

6.1. Solutions for Equivalence Principle 

The solutions presented here use 15,554 ranges from March 1970 to April 2004. The 
ranging stations include the McDonald Observatory, the Observatoire de la Cote 
d'Azur, and the Haleakala Observatory. The ranges of the last decade are fit with a 
2 cm weighted rms residual. Planetary tracking data are used to adjust the orbits 
of the Earth and other planets in joint lunar and planetary fits. The planetary data 
analysis does not include a solution parameter for the equivalence principle. 

Among the solution parameters are GM Eartll _|_M 00n , lunar orbit parameters includ- 
ing semimajor axis, Moon-centered retroreflector coordinates, geocentric ranging 
station coordinates, and lunar tidal displacement Love number hi- For additional 
fit parameters see the modeling discussion in Section [5] An equivalence principle 
violation can be solved for in two ways. The first is a parameter for Mq/Mj with 
a dynamical partial derivative generated from numerical integration. The second 
solves for a coefficient of cosD in the l unar ra nge, a one-term representation. The 
latter approach was used in two paperi^SES Du t, the more sophisticated dynam- 
ical parameter is used in more recent JPL publications, namely Refs. [8l fllll Both 
approaches are exercised here to investigate consistency. 

Five equivalence principle solutions are presented in Table [T] as EP 1 to EP 5. 
Each of these solutions includes a standard set of Newtonian parameters in addition 
to one or more equivalence principle parameters. In addition, the EP solution is 
a comparison case which does not solve for an equivalence principle parameter. 
The solution EP 1 solves for the (Mq/Mi) parameter using the integrated partial 
derivative. That parameter is converted to the coefficient of cos D in radial distance 
using the factor S = -2.9 x 10 13 mm in Eq. ((TTJ from subsection [331 The EP 2 
case solves for coefficients of cos D and sin D in distance using a geometrical partial 
derivative. Solution EP 3 solves for the (Mq/Mj) parameter along with coefficients 
of cos D and sin D in distance. The EP 4 solution constrains the Sun/ (Earth+Moon) 
and Earth/Moon mass ratios. The EP 5 solution uses the mass constraints and also 
constrains the lunar hi- 

The values in Table [2] are corrected for the solar radiation pressure perturbation 
as computed by VokrouhlickyP^ See the modeling subsection !5 . 1 . 9l on solar radiation 
pressure for a further discussion. For the EP 3 case, with two equivalence principle 
parameters, the sum of the two cosD coefficients in Table Q] is —0.6 ± 4.2 mm and 
that sum corrects to 3.1 ±4.2 mm, which may be compared with the four entries in 
Tabled 
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Table 1. Five solutions for the equivalence principle. 



Solution (M G /Mi) 
solution, 
x icr 13 



conversion 
(Ma /Mi) -> coef 
mm 



coef cosD 
solution, 



coef s'mD 
solution, 



bun /(Earth+Moon) 



EP 
EP 1 
EP 2 
EP 3 
EP 4 
EP 5 



0.30 ± 1.42 

0.79 ± 6.09 
0.21 ± 1.30 
-0.11 ± 1.30 



328900.5596 ±0.0011 

-0.9 ±4.1 328900.5595 ± 0.0012 

-0.5 ±4.2 0.9 ±2.1 328900.5596 ± 0.0012 

-2.3 ±17.7 1.7 ±17.8 0.9 ±2.1 328900.5595 ± 0.0012 

-0.6 ± 3.8 328900.5597 ± 0.0007 

0.3 ±3.8 328900.5597 ± 0.0007 



Table 2. Solutions for the equivalence principle corrected for solar radiation pressure. 



Solution 




(M G / Mi) 




(M G /M/) -» coef 


coef cos D 


coef sin D 






solution 




conversion 


solution 


solution 










mm 


mm 


mm 


EP 1 


(- 


0.96 ± 1.42) X 10" 


13 


2.8 ±4.1 






EP 2 










3.1 ±4.2 


0.9 ±2.1 


EP 4 


(- 


1.05 ± 1.30) x 10" 


13 


3.0 ±3.8 






EP 5 


(- 


1.37 ±1.30) x 10" 


13 


4.0 ±3.8 







The equivalence principle solution parameters in Tables [TJand [5] are within their 
uncertainties for all cases except EP 5 in Table [H and that value is just slightly 
larger. Also, the EP 2 coefficient of cos D agrees reasonably for value and uncertainty 
with the conversion of the Mq /Mi parameter of the EP 1 , EP 4 and EP 5 solutions 
to a distance coefficient. For the EP 3 solution, the sum of the converted Mq/Mi 
coefficient and the cos D coefficient agrees with the other solutions in the two tables. 
There is no evidence for a violation of the equivalence principle and solutions with 
different equivalence principle parameters are compatible. 

The difference in uncertainty between the sin D and cos D components of both 
the EP 2 and EP 3 solutions is due to the nonuniform distribution of observations 
with respect to D, as illustrated in Figures [7^, and [8^,. The sinD coefficient is 
well determined from observations near hrst and last quarter Moon, but the cos D 
coefficient is weakened by the decrease of data toward new and full Moon. 

The EP 3 case, solving for Mg/Mi along with cosD and sinD coefficients, is 
instructive. The correlation between the Mg/Mi and cos£> parameters is 0.972 so 
the two quantities are nearly equivalent, as expected. The uncertainty for the two 
equivalence principle parameters increases by a factor of four in the joint solution, 
but the solution is not singular, so there is some ability to distinguish between 
the two formulations. The integrated partial derivative implicitly includes terms at 
frequencies other than the D argument (Nordtvedt, private communication, 1996) 
and it will also have some sensitivity to the equivalence principle influence on lunar 
orbital longitude. The equivalence principle perturbation on lunar orbital longitude 
is about twice the size of the radial component and it depends on sinZ). The ratio 
of Earth radius to lunar semimajor axis is Re/ & ~ 1/60.3, the parallax is about 1°, 
so the longitude component projects into range at the few percent level. 
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The uncertainties in the EP 3 solution can be used to check the theoretical 
computation of the coefficients S, which multiplies A(Mg/Mj), and Co, which 
are associated with the cosD radial perturbation (subsection 3.3). Given the high 
correlation, a first approximation of S = —2.92 x I0 13 mm is given by the ratio of 
uncertainties, and our knowledge that it must be negative. A more sophisticated 
estimate of S = —2.99 x I0 13 mm comes from computing the slope of the axis of the 
uncertainty ellipse for the two parameters. Using expression Eq. ([7]) for the difference 
in self energies of the Earth and Moon, the two preceding values give Ar = 13.0 m 
77 cos D and Ar = 13.3 m 77 cos D, respectively. For comparison, the theoretical 
computations of Ref. 1711 give S — —2.9 x 10 13 mm and Ar = 12.8 m 77 cos D, 
Damour and Vokrouhlicky-S^ give S = —2.9427 x 10 13 mm, corresponding to Ar = 
13.1 m 77 cos D, and Nordtvedt and VokrouhlickjE^ give S — —2.943 x 10 13 mm and 
Ar = 13.1 m 77 cos D. The numerical results here are consistent with the theoretical 
computations within a few percent. 

The EP I solution serves as an example for correlations. The correlation of 
Mq/Mi with both GAf Earth +Moon and osculating semimajor axis (at the 1969 epoch 
of the integration) is 0.46. GM and mean semimajor axis are connected through 
Kepler's third law given that the mean motion is very well determined. The product 
of mean semimajor axis and mean eccentricity is well determined and the correlation 
of Mq I Mi with osculating eccentricity is 0.45. The correlation with the Earth-Moon 
mass ratio is 0.26. 

The value of GMEarth+Moon is important for the equivalence principle solutions. 
The Sun's GM is defined in units of AU 3 /day 2 so GAf Eart h+Moon in those same units 
may be expressed as the mass ratio Sun/(Earth+Moon) as is done in Table [1] The 
Sun/(Earth+Moon) mass ratio is a solution parameter in EP through EP 3. The 
solutions marked EP 4 and EP 5 use a value derived from sources other than LLR. 
The Sun/(Earth+Moon) mass rat io is fixed at a value, with uncertainty, based 
on GM(Earth) from Ries et alP^H and an Earth/Moon mass ratio of 81.300570 ± 
0.000005 from Konopliv et alPSl. The uncertainty for Mq/Mi is improved somewhat 
for solution EP 4. With a fixed GM, the correlation with semimajor axis becomes 
small, as expected, but the correlation with the lunar hi is now 0.42 and the hi 
solution value is 0.044 ± 0.007. For comparison, solution EP 1 had a correlation of 
—0.01 and a solution value of 0.043 ± 0.009. The solution EP 5 adds the lunar Love 
number hi to the cons trained values using hi = 0.0397 from the model calculations 
of Williams et alP^ A realistic model hi uncertainty is about 15%, close to the 
EP 4 solution value, and the Mq/Mi uncertainty is virtually the same as in the EP 
4 solution. All solutions use a model Love number li value constrained to 0.0106. 
Considering the difficulty of precisely comparing uncertainties between analyses 
of different data sets, the gains for the last two constrained equivalence principle 
solutions are modest at best. 

Five solutions presented in this subsection have tested the equivalence principle. 
They do not show evidence for a significant violation of the equivalence principle. 
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6.2. Spectra - Searching for Signatures in the Residuals 

Part of the LLR data analysis is the examination of post-fit residuals including the 
calculation of overall and annual rms, a search for signatures at certain fundamental 
periods, and spectra over a spread of frequencies. Direct examination of residuals 
can reveal some systematic effects but spectra of residuals, appropriately weighted 
for their uncertainties, can expose subtle effects. 

First consider the baseline solution EP without an equivalence principle pa- 
rameter. The distribution of observations vs D has been shown in the histograms 
of Figures [7^ and UK- The last decade of mean weighted residuals vs D is presented 
in Figure [7t> and all of the data is plotted in Figure [8}d. If an equivalence principle 
violation were present it would look like a cosine. No such signature is obvious and 
a fit to the residuals gives a 1 mm amplitude, which is insignificant. 

The LLR data are not evenly spaced or uniformly accurate so aliasing will be 
present in the spectra. Here, a periodogram is computed by sequentially solving for 
sine and cosine components at equally spaced frequencies corresponding to periods 
from 18 years (6585 d) to 6 d. Figure [TUk shows the amplitude spectrum of the 
weighted post-fit residuals for the baseline solution. Nothing is evident above the 
background at the 29.53 day synodic period (frequency #223), which is consistent 
with the results of Table[TJ There are two notable features: a 3.6 mm peak at 1 yr and 
a broad increase at longer periods. There are several uncompensated effects which 
might be contributing at 1 yr including loading effects on the Earth's surface height 
due to seasonal atmosphere and groundwater changes, and "geocenter motion," the 
displacement of the solid body (and core) of the Earth with respect to the overall 
center of mass due to variable effects such as oceans, groundwater and atmosphere. 
Averaged over more than 1000 frequencies the spectrum's background level is 1.2 
mm. Broad increases in the background near 1 month, 1/2 month, and 1/3 month 
etc, are due to aliasing. 

For comparison, an equivalence principle signature was deliberately forced into 
another least-squares solution. A finite A(Ma/Mi) value of 1.5 x 10 -12 , an order of 
magnitude larger than the uncertainty of the EP 1 solution of Tables [T] and [2j was 
constrained in a multiparameter least-squares solution. The standard solution pa- 
rameters were free to minimize the imposed equivalence principle signature as best 
they could. Notably, GAf Earth+M oon and the Earth/Moon mass ratio were distorted 
from normal values by 5 and 3 times their realistic uncertainties, respectively, and 
the correlated orbit parameters also shifted by significant amounts. The overall (35 
year) weighted rms residual increased from 2.9 to 3.1 cm. Figure flOb shows the 
spectrum of the residuals. The two strongest spectral lines, 13 mm and 10 mm, 
are at the D and 3D frequencies, respectively. Detailed examination also shows 
weaker features, a 5D line and mixes of integer multiples of the D frequency with 
the monthly and annual mean anomaly frequencies. The expected equivalence prin- 
ciple signature of 44 mm cos D has been partly compensated by the least-squares 
adjustment of parameters for GM and other quantities. Note that the ratio of the 
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Fig. 10. (a) Spectrum of post-fit residuals without EP solution parameter, (b) Spectrum of resid- 
uals when a A(Mg/M/) value of 1.5 X 10 — 12 is forced into the solution. Frequency #18 corresponds 
to 1 year, #223 is synodic month, and #239 is anomalistic month. 



13 mm peak to the 1.2 mm background is compatible with the ratio of 44 mm (or 
1.5 x 10 -12 ) to the equivalence principle uncertainty of 4.2 mm (or 1.4 x 10~ 13 ) in 
Tables [l]and [2l The spectral amplitudes are computed one frequency at a time, but 
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Table 3. 


Lunar orbit. 




Mean distance 


(r) 


385,000.5 km 


Semimajor axis 


a = l/(l/r> 


384,399.0 km 


Eccentricity 


e 


0.0549 


Inclination to ecliptic plane 


i 


5.145° 



if the amplitudes of cos D and cos 3D are simultaneously fit to the post-fit residuals 
(not to the original ranges) then one gets 34 cos I? + 18 cos 31? in mm. This combi- 
nation would be largest near new Moon, where there are no observations, and near 
full Moon, where there are very few accurate observations. The spectrum for the 
baseline solution in Figure [TUb shows no such lines. In this figure the ~ 3 mm peaks 
near 1 month and 1/3 month are at unassociated periods. 

In summary, a post-fit residual spectrum of baseline solution EP without an 
equivalence principle parameter shows no evidence of any equivalence principle vio- 
lation. Manipulation shows that while a systematic equivalence principle signature 
can be diminished by adjusting other parameters during the least-squares solution, 
that compensation is only partly effective and a systematic effect cannot be elim- 
inated. It is also seen that the parameter uncertainties and correlations from the 
least-squares solutions are in reasonable agreement with the experience based on 
the spectra. 

6.3. Classical Lunar Orbit 

The JPL analyses use numerical integrations for the orbit and dynamical partial 
derivatives. However, Keplerian elements and series expansions for the orbit give 
insight into the solution process. 

The Keplerian elements and mean distance of the Moon are summarized in Table 
131 Note that the inclination is to the ecliptic plane, not the Earth's equator plane. 
The lunar orbit plane precesses along a plane which is close to the ecliptic because 
solar perturbations are much more important than the Earth's J2 perturbation. A 
time average is indicated by (...). 

Various lunar orbital angles and periods are summarized in Table 2J These 
are mean angles represented by smooth polynomials. The solar angles with annual 
periods are V for mean anomaly (the same as the mean anomaly of the Earth-Moon 
center of mass) and V for mean longitude (180° different from the mean longitude 
of the Earth-Moon center of mass). 

The lunar orbit is strongly perturbed by the Sun. Chapront-Touze and Chapront 
have developed an accurate series using computer techniques. From that series (see 
[TBI [TP")) a few large terms for the radial coordinate (in kilometers) are 

r = 385001 - 20905 cos I - 3699 cos(2D - /) - 2956 cos 2D - 

- 570cos2Z + 246cos(2Z-2.D) + ... + 109cosL> + ... (19) 

The constant first term on the right-hand side is the mean distance (somewhat larger 
than the semimajor axis), the I and 21 terms are elliptical terms, and the remaining 
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Table 4. Lunar angles. 



Angle 


Symbol 


Period 


Mean Longitude 


L 


27.322 d 


Mean Anomaly 


I 


27.555 d 


Mean Argument of Latitude 


F 


27.212 d 


Mean Elongation of Moon from Sun 


D 


29.531 d 


Mean Node 


n 


18.61 yr 


Mean Longitude of Perigee 




8.85 yr 


Mean Argument of Perigee 


Ul 


6.00 yr 



terms are from solar perturbations. The amplitudes of the solar perturbation terms 
depend on the masses of the Earth, Moon, and Sun, as well as the lunar orbit and 
the Earth-Moon orbit about the Sun. The periods of the periodic terms in the order 
given in Eq. (O are 27.555 d, 31.812 d, 14.765 d, 13.777 d, 205.9 d, and 29.531 d, 
so the different terms are well separated in frequency. 

If the equivalence principle is violated, there is a dipole term in the expansion 
of the solar perturbation which gives the cosD term of subsection 13. 3[ see Refs. [67j 
[7T1 When the equivalence principle is satisfied the dipole term has zero coefficient. 
There is a classical cosl? term which arises from the octupole (P3) term in the 
expansion and that gives the 109 km amplitude in the series expansion for orbital 
r, Eq. CEU). 

The JPL Lunar Laser Ranging analyses use numerically integrated orbits, not 
series expansions (see [20] for the polynomial expressions for lunar angles and an 
LLR data analysis with a higher reliance on analytical series). The uncertainty of 
the solar perturbation corresponding to the classical cos D term is very small and 
is included in the final Mq/Mj and amplitude uncertainties of the EP 1, EP 2, and 
EP 3 solutions of Tables Q] and [21 since mass and orbit quantities are also solution 
parameters in those least-squares solutions. 

6.4. Separation of the Equivalence Principle Signature 

The equivalence principle solution parameter, whether Mq/Mi or cosD, is signif- 
icantly correlated with GM of the Earth-Moon system and lunar semimajor axis. 
The mean motion of the Moon is very well determined from the observations so Ke- 
pler's third law strongly relates the GM and mean semimajor axis. The correlation 
between GM and cos D is related to the uneven distribution of observations for the 
angle D (Figures [7^,00- The relation between the equivalence principle, GM and 
the D distribution has been extensively discussed by Nordtvedfr^. Some additional 
effects are briefly described by Anderson and Williams^. This subsection discusses 
the consequence of the D distribution and other effects. 

The range may be derived from the vector Eq. (|14[) . The scalar range equation 
may be approximated as 

p « r - (r • R stn ) + (f • Rrfx) - (R stn • RrfO/r + ... (20) 

The extended series expansion of the range equation is complicated, but with some 
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consideration a few terms may be selected which are relevant to the equivalence 
principle solution. The dot product between the orbital radius and the station vec- 
tor involves large near daily and monthly variations and in solutions the station 
coordinates separate well from the other parameters. Because of the good separa- 
tion, this dot product will not be considered further here. The series for orbital 
radius r is given in Eq. (fl"9l) . The mean distance for an elliptical orbit is given by 
o(l + e 2 /2), where for the Moon a = 384,399 km and e = 0.0549 (see Table©. 
The terms with mean anomaly I in the arguments have coefficients that depend 
on eccentricity. The coefficient of the 2D term is scaled by the semimajor axis and 
while the coefficient has sensitivity to other parameters the semimajor axis scaling 
is the primary concern here. The mean anomaly dependent terms have periods quite 
different from D and will not be considered further. Considering the dot product 
between reflector and orbit radius, the term of interest is Xu\, where X is the reflec- 
tor component toward the mean Earth direction, expressed in the body-referenced 
frame, and the expansion for the x componen t of t he unit vector from Moon center 
to Earth center in the same frame is 

mi w 0.99342 + 0.00337 cos 2F + 0.00298 cos 21 + 

+ 0.00131 cos2L> - 0.00124cos(2; - 2D) + ... (21) 

The angle F is the polynomial for the mean argument of latitude, the lunar angle 
measured from the node of the orbit on the ecliptic plane. This angle is associated 
with the tilts of the orbit plane and the lunar equator plane to the ecliptic plane 
and it has a period of 27.212 days (Tabled}. 

With the above considerations the relevant combination of terms is 

NcosD + o(1.00157- 0.00769 cos 2D) - Xu x - (R stn • R rf i)/r, (22) 

where the first term represents an equivalence principle violation and N, a, and X 
are to be determined from the data. The linear combination 1.0016a — 0.9934X is 
better determined by two orders-of-magnitude than either a or X. The separation of 
the different solution parameters is aided by the time variation of their multiplying 
functions in Eq. (|22|) . The periodic 2D term provides one way to separate X and 
a. If the angle D were uniformly distributed, then the D and 2D terms would 
be distinct. The nonuniform distribution of D (Figures [7K and |HK) weakens the 
separation of the two periodicities and causes N, a [and GM Eartll _|_M 00n ] , and X to 
be correlated. The separation of X is aided by the periodic terms in Ui, such as 
the two half month terms with arguments 2F and 21, as well as the dot product 
between the station and reflector vectors, where R s tn/a = 1/60.3 sets the scale for 
daily and longer period terms. 

A good equivalence principle test is aided by a) a good distribution of angle 
D, b) a good distribution of orbit angles I and F, which is equivalent to a good 
distribution of orientations of the Moon's x axis with respect to the direction to the 
Earth (optical librations), and c) a wide distribution of hour angles and declinations 
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of the Moon as seen from the Earth. Of these three, the first is the hardest for LLR 
to achieve for the reasons discussed in subsection 14.21 on selection effects. 

7. Derived Effects 

The solution EP 1 matches the EP test published in Ref. 11111 The data set of this 
paper has only one data point more than the data set of the published case. Several 
consequences can be derived from the equivalence principle test including a test of 
the strong equivalence principle and PPN parameter 0. 

7.1. Gravity Shielding - the Majorana Effect 

The possibility that matter can shield gravity is not predicted by modern theories 
of gravity, but it is a recurrent idea and it would cause a violation of the equivalence 
principle test. Consequently, a brief discussion is given in this subsection. 

The idea of gravity shielding goes back at least as far as to the original paper by 
MajoranaP^He proposed that the inverse square law of attraction should include an 
exponential factor exp(— h J p(s)ds) which depends on the amount of mass between 
attracting mass elements and a universal constant h. If mass shields gravity, then 
large bodies such as the Moon and Earth will partly shield their own gravitational 
attraction. The observable ratio of gravitational mass to inertial mass would not 
be independent of mass, which would violate the equivalence principle. Russell^ 
realized that the large masses of the Earth, Moon and planets made the observations 
of the orbits of these bodies and the asteroid Eros a good test of such a possibility. 
He made a rough estimate that the equivalence principle was satisfied to a few 
parts per million, which was much smaller than a numerical prediction based on 
Majorana's estimate for h. 

Majorana gave a closed form expression for a sphere's gravitational to inertial 
mass ratio. For weak shielding a simpler expression is given by the linear expansion 
of the exponential term 



where / is a numerical factor, p is the mean density, and R is the sphere's radius. 
For a homogeneous sphere Majorana and Russell give / = 3/4. For a radial density 
distribution of the form p(r) = p(0)(l — r 2 / R 2 ) n Russell - derives / = (2n + 
3)7(12n+ 12). 



EckhardtF^ used an LLR test of the equivalence principle to set a modern limit 
on gravity shielding. That result is updated as follows. The uniform density approx- 
imation is sufficient for the Moon and fRp — 4.4 x 10 8 gm/cm 2 . For the Earth we 
use n ss 0.8 with Russell's expression to get fRp = 3.4 x 10 9 gm/cm 2 . Using the 
difference —3.0 x 10 9 gm/cm 2 h along with the LLR EP 1 solution from Table [2] 
for the difference in gravitational to inertial mass ratios gives h = (3 ± 5) x 10~ 23 



hfRp, 

Mi 



(23) 
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cm 2 /gm. The value is not significant compared to the uncertainty. To give a sense 
of scale to the uncertainty, for the gravitational attraction to be diminished by 1/2 
would require a column of matter with the density of water stretching at least half 
way from the solar system to the center of the galaxy. The LLR equivalence principle 
tests give no evidence that mass shields gravity and the limits are very strong. 

7.2. The Strong Equivalence Principle 

The total equivalence principle results for the Earth-Moon system have been given 
in Tabled This test is a strong result in its own right. The total equivalence principle 
is the sum of contributions from the WEP, which depends on composition, and the 
SEP, which depends on gravitational self energy. This subsection extracts a result 
for the SEP by using WEP results from laboratory experiments at the University 
of Washington. 

Experiments by several groups have tested the WEP. Several of these experi- 
ments with different test body compositions were compared in order to limit the 
WEP effect on the Earth-Moon pair to 10~ 12 , see Refs. 121 111 Recen t la boratory 
investigations have synthesized the composition of the Earth and MoorJ^El by using 
test bodies which simulate the composition of core and mantle materials. These 
WEP results are an order-of-magnitude more accurate. 

The most abundant element in the Earth is oxygen, followed by iron (30 weight 
%), silicon and magnesiumlSSl For the Moon, iron is in fourth place with about 1/3 
of the Earth's abundance. The composition of the mantles of the Earth and Moon 
are similar, though there are differences (e.g. the Moon lacks the lower tempera- 
ture volatiles such as water). Iron and nickel are the heaviest elements which are 
abundant in both bodies. Hence the difference in iron abundance, and associated 
siderophile elements, between the Earth and Moon is the compositional difference 
of most interest for the WEP. 

The Earth has a massive core (~l/3 by mass) with iron its major constituent and 
nickel and sulfur lesser components. Several lines of evidence indicate that the Moon 
has a small core which is < 2% of its mass: mom ent of inertiePS, induced magnetic 
dipole momenljSSI, and rotational dynamic d^*! The lunar core is presumed to be 
dominated by iron, probably alloyed with nickel and possibly sulfur, but the amount 
of information on the core is modest and evidence for composition is indirect. In 
any case, most of the Fe in the Moon is in minerals in the thick mantle while for 
the Earth most of the Fe is in the metallic core. For an example of lunar models see 

Ref . EU E2 

For consideration of the WEP the iron content is the important difference in 
composition between the Earth and Moon. Among the elements present at > 1 
weight %, iron (and nickel for the Earth) have the largest atomic weights and 
numbers. The two University of Washington test bodies reproduce the mean atomic 
weights and mean number of neutrons for the core material of the Earth (and 
probably the Moon) and both bodies' mantles. 
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The Baessler et alPand Adelberger^ analyses use 38.2 % for the fraction of mass 
of Fe/Ni core material in the whole Earth and 10.1 % for the fraction in the Moon. 
The difference in the experimental accelerations of the two test bodies is converted 
to the equivalent (WEP) difference in the acceleration of the Earth and Moon by 
multiplying by the difference (0.281). Since the iron contents of the Earth and 
Moon are uncertain by a few percent, the effect of composition uncertainties is an 
order-of-magnitude less than the derived acceleration difference. The AdelbergeJ^ 
result for the relative acceleration is given as (1.0 ± 1.4 ± 0.2) x 10 -13 , where the 
first uncertainty is for random errors and the second is for systematic errors. We 
combine the systematic and random uncertainties and use 



(Mg\ _ (Ma 



Mr 



M. 



(1.0 ± 1.4) x 10 



-13 



(24) 



The strong equivalence principle test comes from combining solution EP 1 of Tabled] 
with the above WEP result. 



= (-2.0 ±2.0) x 10~ 13 . (25) 

SEP 

This combination of the LLR determination of the equivalence principle and the 
laboratory test of the weak equivalence principle provides the tightest constraint on 
the strong equivalence principle. 




7.3. PPN Beta 

The test for a possible violation of the strong equivalence principle, the equivalence 
principle due to self-energy, is sensitive to a linear combination of PPN parameters. 
For conservative theories this linear relation is r\ = 4/3—7—3, given by Eq. ((3|). Using 
a good experimental determination of PPN 7, the SEP result can be converted into 
a result for PPN ft. 

The test for any violation of the strong equivalence principle is sensitive to a 
linear combination of PPN quantities. Considering only PPN ft and 7, divide the 
SEP determination of Eq. (f25f by the numerical value from Eq. ([7|) to obtain 

7^ = 4/3-7-3= (4.4 ±4.5) x 10~ 4 . (26) 

This expression would be null for general relativity, hence the small value is consis- 
tent with Einstein's theory. 

The SEP relates to the non-linearity of gravity (how gravity affects itself) , with 
the PPN parameter ft representing the degree of non-linearity. LLR provides great 
sensitivity to ft, as suggested by the strong dependence of 77 on ft in Eqs. ([3]) and 

(23). 

An acc ura te result for 7 has been determined by the Cassini spacecraft 
experiment Using high-accuracy Doppler measurements, the gravitational time 
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delay allowed 7 to be determined to the very high accuracy of 7 — 1 = (2.1 ± 2.3) x 
10~ 5 . This value of 7, in combination with 77, leads to a significant improvement in 
the parameter (3: 

p- 1 = (1.2 ± 1.1) x nr 4 . (27) 

We do not consider this result to be a significant deviation of (3 from unity. 

The PPN parameter (3 has been determined by combining the LLR test of the 
equivalence principle, the laboratory results on the WEP, and the Cassini spacecraft 
determination of 7. The uncertainty in (3 is a dramatic improvement over earlier 
results. The data set for the solutions in this chapter differs by only one point from 
that used in Ref. 1 1 1 1 1 Consequently, the equivalence principle solution EP 1, and 
the derived result above for the strong equivalence principle, 77 and (3 are virtually 
the same as for the publication. 

8. Emerging Opportunities 

It is essential that the acquisition of new LLR data continue in the future. Cen- 
timeter level accuracies are now achieved, and a further improvement is expected. 
Analyzing improved data would allow a correspondingly more precise determination 
of gravitational physics and other parameters of interest. In addition to the existing 
LLR capabilities, there are two near term possibilities that include the construc- 
tion of the new LLR stations and development and deployment of cither new sets 
of passive laser cornercube retroreflectors or active laser transponders pointed at 
Earth or both of these instruments. 

In this Section we will discuss both of these emerging opportunities - the new 
LLR station in New Mexico and new LLR instruments on the Moon - for near term 
advancements in gravitational research in the solar system. 

8.1. New LLR Data and the APOLLO facility 

LLR has remained a viable experiment with fresh results over 35 years because 
the data accuracies have improved by an order of magnitude (see Figure [5J . A 
new LLR station should provide another order of magnitude improvement. The 
Apache Point Observatory Lunar Laser-ranging Operation (APOLLO) is a new 
LLR effort designed to achieve millimeter range precision and corresponding order- 
of-magnitude gains in measurements of fundamental physics parameters. Using a 3.5 
m telescope the APOLLO facility will push LLR into the regime of stronger photon 
returns with each pulse, enabling millimeter range precision to be achievedP^Hini 
An advantage that APOLLO has over current LLR operations is a 3.5 m astro- 
nomical quality telescope at a good site. The site in southern New Mexico offers 
high altitude (2780 m) and very good atmospheric "seeing" and image quality, with 
a median image resolution of 1.1 arcseconds. Both the image sharpness and large 
aperture combine to deliver more photons onto the lunar retroreflector and receive 
more of the photons returning from the reflectors, respectively. Compared to current 
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operations that receive, on average, fewer than 0.01 photons per pulse, APOLLO 
should be well into the multi-photon regime, with perhaps 1-10 return photons 
per pulse, depending on seeing. With this signal rate, APOLLO will be efficient at 
finding and tracking the lunar signal, yielding hundreds of times more photons in 
an observation than current operations deliver. In addition to the significant reduc- 
tion in random error (l/y/~N reduction), the high signal rate will allow assessment 
and elimination of systematic errors in a way not currently possible. This station 
is designed to deliver lunar range data accurate to one millimeter. The APOLLO 
instrument started producing useful ranges in 2006 , thereby, ini tiating the regular 
delivery of LLR data with much improved accuracy! I I I I 

The high accuracy LLR station installed at Apache Point should provide major 
opportunities (see Refs. I6"2l 11101 IB"4l for details). The APOLLO project will push 
LLR into the regime of millimctric range precision which translates into an order- 
of- magnitude improvement in the determination of fundamental physics parameters. 
An Apache Point 1 mm range accuracy corresponds to 3 x 10 -12 of the Earth-Moon 
distance. The resulting LLR tests of gravitational physics would improve by an 
order of magnitude: the Equivalence Principle would give uncertainties approaching 
10~ 14 , tests of general relativity effects would be <0.1%, and estimates of the relative 
change in the gravitational constant would be 0.1% of the inverse age of the universe. 
This last number is impressive considering that the expansion rate of the universe 
is approximately one part in 10 10 per year. Therefore, the gain in our ability to 
conduct even more precise tests of fundamental physics is enormous, thus this new 
instrument stimulates development of better and more accurate models for the LLR 
data analysis at a mm-level. 



8.2. New retroreflectors and laser transponders on the Moon 

There are two critical factors that control the progress in the LLR-enabled science 
- the distribution of retroreflectors on the lunar surface and their passive nature. 
Thus, the four existing arrays^l are distributed from the equator to mid-northern 
latitudes of the Moon and are placed with modest mutual separations relative to 
the lunar diameter. Such a distribution is not optimal; it limits the sensitivity of the 
ongoing LLR science investigations. The passive nature of reflectors causes signal 
attenuation proportional to the inverse 4th power of the distance traveled by a laser 
pulse. The weak return signals drive the difficulty of the observational task; thus, 
only a handful of terrestrial SLR stations are capable of also carrying out the lunar 
measurements, currently possible at cm-level. 

The intent to return to the Moon was announced in January 2004. NASA is 
planning to return to the Moon in 2009 with Lunar Reconnaissance Orbiter, and 
later with robotic landers, and then with astronauts in the next decade. The return 
to the Moon provides an excellent opportunity for LLR, particularly if additional 
LLR instruments will be placed on the lunar surface at more widely separated 
locations. Due to their potential for new science investigations, these instruments 
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are well justified. 

8.2.1. New retroreflector arrays 

Future ranging devices on the Moon might take two forms, namely passive retrore- 
flectors and active transponders. The advantages of passive retroreflector arrays are 
their long life and simplicity. The disadvantages are the weak returned signal and 
the spread of the reflected pulse arising from lunar librations, which can change 
the retroreflector orientation up to 10 degrees with respect to the direction to the 
Earth. 

Range accuracy, data span, and distributions of earth stations and retroreflectors 
are important considerations for future LLR data analysis. Improved range accuracy 
helps all solution parameters. Data span is more important for some parameters, 
e.g. change in G, precession and station motion, than others. New retroreflectors 
optimized for pulse spread, signal strength, and thermal effects, will be valuable at 
any location on the moon. 

Overall, the separation of lunar 3-dimensional rotation, the rotation angle and 
orientation of the rotation axis (also called physical librations), and tidal displace- 
ments depends on a good geographical spread of retroreflector array positions. The 
current three Apollo sites plus the infrequently observed Lunokhod 2 are close to 
the minimum configuration for separation of rotation and tides, so that unexpected 
effects might go unrecognized. A wider spread of retroreflectors could improve the 
sensitivity to rotation/orientation angles and the dependent lunar science param- 
eters by factors of up to 2.6 for longitude and up to 4 for pole orientation. The 
present configuration of retroreflector array locations is quite poor for measuring 
lunar tidal displacements. Tidal measurements would be very much improved by 
a retroreflector array near the center of the disk, longitude and latitude 0, plus 
arrays further from the center than the Apollo sites. 

Lunar retroreflectors are the most basic instruments, for which no power is 
needed. Deployment of new retroreflector arrays is very simple: deliver, unfold, point 
toward the Earth and walk away. Retroreflectors should be placed far enough away 
from astronaut /moonbase activity that they will not get contaminated by dust. 
One can think about the contribution of smaller retroreflector arrays for use on 
automated spacecraft and larger ones for manned missions. One could also benefit 
from co-locating passive arrays and active transponders and use a few LLR capable 
stations ranging retroreflectors to calibrate the delay vs. temperature response of 
the transponders (with their more widely observable strong signal). 

8.2.2. Opportunity for laser transponders 

LLR is one of the most modern and exotic observational disciplines within astrom- 
ctry, being used routinely for a host of fundamental astronomical and astrophysical 
studies. However, even after more than 30 years of routine observational operation, 
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LLR remains a non-trivial, sophisticated, highly technical, and remarkably chal- 
lenging task. Signal loss, proportional to the inverse 4th power of the Earth-Moon 
distance, but also the result of optical and electronic inefficiencies in equipment, ar- 
ray orientation, and heating, still requires that one observe mostly single photoelec- 
tron events. Raw timing precision is some tens of picoseconds with the out-and-back 
range accuracy being approximately an order of magnitude larger. Presently, we are 
down to sub-cm lunar ranging accuracies. In this day of routine SLR operations, it 
is a sobering fact to realize that ranging to the Moon is many orders of magnitude 
harder than to an Earth-orbiting spacecraft. Laser transponders may help to solve 
this problem. Simple time-of-flight laser transponders offer a unique opportunity 
to overcome the problems above. Although there are great opportunities for sci- 
entific advances provided by these instruments, there are also design challenges as 
transponders require power, precise pointing, and thermal stability. 

Active laser transponders on the lunar surface are attractive because of the 
strong return and insensitivity to lunar orientation effects. A strong return would 
allow artificial satellite ranging stations to range the Moon. However, transponders 
require development: optical transponders detect a laser pulse and fire a return pulse 
back toward the Earth.ES They give a much brighter return signal accessible to 
more stations on Earth. Active transponders would require power and would have 
more limited lifetimes than passive reflectors. Transponders might have internal 
electronic delays that would need to be calibrated or estimated, since if these delays 
were temperature sensitive that would correlate with t he SE P test. Transponders 
can also be used to good effect in asynchronous wherein the received 

pulse train is not related to the transmitted pulse train, but the transponder unit 
records the temporal offsets between the two signals. The LLR experience can help 
determine the optimal location on the Moon for these devices. 

In addition to their strong return signals and insensitivity to lunar orientation 
effects, laser transponders are also attractive due to their potential to become in- 
creasingly important part of space exploration efforts. Laser transponders on the 
Moon can be a prototype demonstration for later laser ranging to Mars and other ce- 
lestial bodies to give strong science returns in the areas similar to those investigated 
with LLR. A lunar installation would provide a valuable operational experience. 



9. Summary 

In this paper we considered the LLR tests of the equivalence principle (EP) per- 
formed with the Earth and Moon. If the ratio of gravitational mass to inertial mass 
is not constant, then there would be profound consequences for gravitation. Such 
a violation of the EP would affect how bodies move under the influence of gravity. 
The EP is not violated for Einstein's general theory of relativity, but violations are 
expected for many alternative theories of gravitation. Consequently, tests of the EP 
are important to the search for a new theory of gravity. 

We considered the EP in its two forms (Sec. [3]); the weak equivalence princi- 
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pie (WEP) is sensitive to composition while the strong equivalence principle (SEP) 
considers possible sensitivity to the gravitational energy of a body. The main sensi- 
tivity of the lunar orbit to the equivalence principle comes from the acceleration of 
the Earth and Moon by the Sun. Any difference in those accelerations due to a fail- 
ure of the equivalence principle causes an anomalous term in the lunar range with 
the 29.53 d synodic period. The amplitude would be proportional to the difference 
in the gravitational to inertial mass ratios for Earth and Moon. Thus, lunar laser 
ranging is sensitive to a failure of the equivalence principle due to either the WEP 
or the SEP. In the case of the SEP, any violation of the equivalence principle can 
be related to a linear combination of the parametrized post-Newtonian parameters 
(3 and 7. 

We also discussed the data and observational influences on its distribution 
(Sec. |4j. The evolution of the data from decimeter to centimeter quality fits is 
illustrated. The LLR data set shows a variety of selection effects which influence 
the data distribution. Important influences include phase of the Moon, season, dis- 
tance, time of day, elevation in the sky, and declination. For the LLR-enabled EP 
tests, selection with phase of the Moon is an important factor. 

An accurate model and analysis effort is needed to exploit the lunar laser range 
data to its full capability. The model is the basis for the computer code that pro- 
cesses the range data (Sec.[5|). Further modeling efforts will be necessary to process 
range data of millimeter quality. Two small effects for future modeling, thermal 
expansion and solar radiation pressure, are briefly discussed. 

Solutions for any EP violation are given in Section [5] Several approaches to the 
solutions are used as checks. The EP solution parameter can be either a ratio of 
gravitational to inertial masses or as a coefficient of a synodic term in the range 
equation. The results are compatible in value and uncertainty. Because GM^ rt)1+Koon 
correlates with the EP due to lunar phase selection effects, solutions are also made 
with this quantity fixed to a value based on non-LLR determinations of GMEarth 
and Earth/Moon mass ratio. In all, five EP solutions are presented in Tableland 
four are carried forward into Table[2j As a final check, spectra of the post-fit residuals 
from a solution without any EP solution parameter are examined for evidence of any 
violation of the EP. No such signature is evident. The analysis of the LLR data does 
not show significant evidence for a violation of the EP compared to its uncertainty. 
The final result for [(Mg/M^e - {M G /Mi) M ] E p is (-1.0 ± 1.4) x 10~ 13 . 

To gain insight into the lunar orbit and the solution for the EP, short trigono- 
metric series expansions are given for the lunar orbit and orientation which are 
appropriate for a range expansion. This is used to show how the data selection with 
lunar phase correlates the EP solution parameter with GM Earth+M oon- To separate 
these and other relevant parameters, one wishes a good distribution of observations 
with lunar phase, orbital mean anomaly and argument of latitude, and, as seen from 
Earth, hour angle and declination. 

The result for the SEP is derived (subsection 17. 2p from the total value deter- 
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mined by LLR by subtracting the laboratory result for the WEP determined at 
the University of Washington. The Moon has a small core while the Earth has a 
large iron rich core. Both have silicate mantles. The WEP sensitivity of the Moon 
depends most strongly on the difference in iron content between the two bodies. 
The SEP result is [{M G /Mi) E - (M G /Mj) M ] SEP = (-2.0 ± 2.0) x 10~ 13 , which we 
do not consider to be a significant difference from the zero of general relativity. 

The SEP test can be related to the parametrized post-Newtonian (PPN) pa- 
rameters /3 and 7 (subsection 17.3ft . For conservative theories of relativity, one gets 
4/3 — 7 — 3 = (4.4 ± 4.5) x 10~ 4 . The Cassini spacecraft result for 7 allows a value 
for (3 to be extracted. That result is (3 — 1 = (1.2 ± 1.1) x 10~ 4 , which is the most 
accurate determination to date. Again, we do not consider this (3 value to be a 
significant deviation from the unity of general relativity. 

Finally, we discussed the efforts that are underway to extend the accuracies 
to millimeter levels (Sec. [8]). The expected improvement in the accuracy of LLR 
tests of gravitational physics expected with extended data set with existing stations 
and also with a new APOLLO instrument will bring significant new insights to 
our understanding of the fundamental physics laws that govern the evolution of our 
universe. The scientific results are very significant which justifies the nearly 40 years 
of history of LLR research and technology development. 

The lunar laser ranging results in this paper for the equivalence principle, strong 
equivalence principle, and PPN j3 are consistent with the expectations of Einstein's 
general theory of relativity. It is remarkable that general relativity has survived a 
century of testing and that the equivalence principle is intact after four centuries of 
scrutiny. Each new significant improvement in accuracy is unknown territory and 
that is reason for future tests of the equivalence principle. 
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